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The probable mechanisms underlying the development of boundary layer roll circulations 
are studied using wind and temperature profiles measured by the NCAR Electra during the 
stratocumulus phase of the FIRE experiment The expected, or preferred, roll orientations, 
horizontal wavelengths and propagation periods are determined by finding the minimum values 
of the dynamic and thermodynamic forcing parameters, which here are the eddy Reynolds 
number Re and moist Rayleigh number Ra^ These minimum values depend on the height zj of 

the capping temperature inversion and on the values of the Fourier coefficients of the 
background height-dependent vector wind profile. As input to our nonlinear spectral model, 
descent and ascent runs by the Electra provide good initial estimates of the inversion height and 
the wind profiles. In the first phase of the investigation presented here, a mechanism is said to 
be a probable contributor to the development of roll circulations within die stratocumul us-topped 
boundary layer if the modeled roll orientations and wavelengths agree with their observed 
values. 


Possible mechanisms for roll development are the well-known thermal and inflection 
point instability mechanisms (e.g. Brown, 15)80) and the recently identified shear instability 
mechanism of Haack-Hirschberg (1988). The first mechanism leads to thermal modes 
dominated by one wavenumber in the vertical, while the second two lead to dynamic modes 
requiring two wavenumbers in the vertical. Several low-order spectral models that are capable 
of accepting observed data for determination of the expected modes have been developed at 
Penn State; these include the pure thermal model of Snirer (1986), the pure inflection point 
model of Stensrud and Shirer (1988), and the mixed thermal/dynamic model of 
Haack-Hirschberg (1988). Linear analyses of the solutions to these models produce estimates of 
the preferred values of the orientation angles, horizontal wavelengths, and dimensionless 
frequencies for the above three modes, and these values can be readily compared with 
observations. 

Preliminary results using the 14-coefficient model of Haack-Hirschberg (1988) are 
discussed for the 7 July 1987 Electra Mission 188-A (Flight 5). This mission was flown across a 
sharp cloud boundary that was within a Landsat/SPOT scene (Kloesel et al., 1988). The 
stratocumulus deck was relatively solid in the eastern part of the scene, while there was a rapid 
decrease in cloud cover to scattered cumulus clouds aligned in streets to the west. These cloud 
streets were oriented nearly parallel to the mean wind direction in the layer, which was 

approximately 340°. The hypothesis that roll circulations occurred in both the relatively clear 
and the cloudy regions is investigated using as model input a descent profile obtained in the 
relatively clear air and an ascent profile obtained in the cloudy air. Initial results for the clear air 
case are that the pure inflection point mode is not possible and the pure thermal mode was 

oriented 35° to the right of the mean wind direction. The origin of this unacceptably large 
discrepancy between the observed and modeled results will be investigated further and the 
conclusions reported at the next FIRE workshop. 

The results for the cloudy case are more promising. In Fig.l are shown the measured 
cross-roll U*j and along-roll V^j wind components (jagged curves), together with the fit to the 
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data implied by the model representation (smooth curves). Clearly the limited resolution of the 
14-coefficient model is able to represent the observations rather well. As shown by the solid and 
dashed lines in Fig. 2, two thermal modes are possible in this case. The quadrilateral in Fig. 2 
encloses the estimated range of observed values of Re and Rag}, whore Re ■ IV(zt)Izx/(xv), 

Rag, 3 gzf A^TAvxToow 4 ), IV(z T )l is the wind at the inversion base zj, v and k are the eddy 

viscosity and thermometric conductivity respectively, and AjjT is a combined measure of the sea 

surface/air temperature difference and the environmental lapse rate of equivalent potential 
temperature. The mode whose transition curve passes to the left of the estimated range of 

observed values for Ra,n and Re is associated with orientations that are 10° to the right of the 

mean wind direction (335°); this prefared orientation angle 0 p is 75° south of east, or -75°. The 
expected aspect ratio ap of 0.6 corresponds to a horizontal roll wavelength L * 2z-pfap of 

approximately 2500 m. This value of L is consistent with the wavelengths found from 
preliminary vertical velocity spectra for the horizontal flight leg between the descent and ascent 
profiles. Finally the dimensionless frequency 0)p of 31.7 corresponds to a propagation period of 

approximately 30 minutes. Thus, the model results suggest that a thermally driven roll 
circulation supported the stratocumulus clouds observed on 7 July 1987. 
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Figure 1. Cross-roll U M and along-roll V M wind profiles (jagged curves) measured by the 
NCAR Electra in the stratocumulus region on 7 July 1987 (Flight 5) between 2154:43 
and 2156:49 UTC. The cross-roll, or x-, direction is 65° (~ENE), so that negative values 

correspond to winds having directions less than 335°. The representation of the data by 
the model is shown by the smooth curves. 
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Figure 2. Results from the model analysis. The solid and dashed lines denote the transition 
curves for the two thermal modes supported by the observed wind profile given in Fig. 1. 
The numbers next to these curves denote the preferred orientation angle Pp. the aspect 
ratio 3 p and the dimensionless frequency U) p (in parentheses). The quadrilateral encloses 
the values of the moist Rayleigh number Ra m and Reynolds number Re that are estimated 

to characterize the observed stratocumulus regime. These vertices are given 
approximately by (Ra m ,Re) = (0,70); (10,70); (145,285); (5,285). Here z T - 750 m and 
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v and K vary between 10 m /s and 50 m /s. 
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